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A series of density functional theory (DFT) simulations was performed to investigate the approach, 
contact, and subsequent separation of two atomically flat surfaces consisting of different materials. 
Aluminum (Al) and titanium nitride (TiN) slabs were chosen as a model system representing a 
metal-ceramic interface and the interaction between soft and hard materials. The approach and 
separation were simulated by moving one slab in discrete steps normal to the surfaces allowing 
for electronic and atomic relaxations after each step. Various configurations were analyzed by 
considering (001), (011), and (111) surfaces as well as several lateral arrangements of these surfaces 
at the interface. Several tests were conducted on the computational setup, for example, by changing 
the system size or using different approximations for the exchange correlation functional. The 
performed simulations revealed the influences of these aspects on adhesion, equilibrium distance, 
and material transfer. These interfacial properties depend sensitively on the chosen configuration 
due to distinct bond situations. Material transfer, in particular, was observed if the absolute value 
of the adhesion energy for a given configuration is larger than the energy cost to remove surface 
layers. This result was found to be independent of the employed exchange correlation functional. 
Furthermore, it was shown that a simple comparison of the surface energies of the slabs is not 
sufficient to predict the occurrence of material transfer. 
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I. INTRODUCTION 

Contacts of surfaces at the atomic length scale are 
crucial in many modern applications, from experimental 
techniques such as nanoindentatiorP^ or atomic/friction 
force microscopy (AFM/FFM)PM^ to nanotechnologies 
applied, for example, in nano-/microelectromechanical- 
systems (NEMS/MEMS)Pni The reliability, perfor¬ 
mance, and lifetime of such systems, for example, depend 
sensitively on the interactions between contacting mate¬ 
rials. Furthermore, detailed insights into such contacts 
are of fundamental interest for better comprehension of 
tribological processes, such as nanoscale wear,E^K^ for 
which there is still a lack of understanding due to its 
highly complex nature.!^ 

Metal-ceramic interface^^ are of fundamental and 
technological interest because they exhibit advantages 
of both types of materials such as valuable mechanical 
properti^ high thermal stability, and degradation re¬ 
sistance.^ Hence, such interfaces are important in nu¬ 
merous applications such as communication devices and 
nanoelectronics.l^In this paper the interface between the 
metal Al and the transit ion-metal nitride TiN is inves¬ 
tigated. This interface consists of a soft material and a 
hard material, which simplifies wear processes because 
the softer material is primarily affected. 

Since the 1980s classical molecular dynamics (MD) 
simulations have commonly been applied to nanotribo- 
logical problems (see, e.g., Refs. [24143^ and still consti¬ 
tute a standard tool in numerical atomistic simulations. 
Nevertheless, during the last decade density functional 
theory (DFT) calculations have been increasingly used 


in this field (see, e.g., Refs. [35H46|) and should be seen as 
an extension to the more common computational tools 
in tribology. DFT allows for parameter-free calculations 
and an accurate description of quantum-mechanical sys¬ 
tems and does not depend on empirical potentials. How¬ 
ever, due to computational challenges DFT calculations 
are currently limited to rather small systems of typi¬ 
cally a few hundred atoms. Since DFT has proven to 
yield reliable results for this class of systems 
also employed in this study to analyze the electronic and 
atomic structure of the Al/TiN interfaces and to deter¬ 
mine properties such as adhesion energies. Results ob¬ 
tained with DFT, such as potential-energy curves, can be 
used as an i nput for, e.g., large-scale classical MD simu¬ 
lations Furthermore, the combination of approaches 
such as DFT and MD as well as the continuously increas¬ 
ing available computer power and advances in software 
tools promise the possibility to investigate even larger 
and more realistic systems in the near future. 

Al/TiN and similar interfaces have already bee n in- 
vestigated by various researchers with experimentaP^^^ 
and theoreticaP^M^ methods. Here, however, the em¬ 
phasis lies on a realistic way to simulate the separation 
of the interfaces as well as on a comprehensive discussion 
of interfaces between Al and TiN low-index surfaces. To 
assess this problem, the effects of various configurations 
at the interface as well as approach and subsequent sep¬ 
aration of Al and TiN slabs are analyzed. Various tests 
on, e.g., the effect of adjusted lattice parameters, the 
simulation cell size, and various approximations for the 
exchange correlation functional in DFT are carried out. 
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II. COMPUTATIONAL DETAILS 
A. Density Functional Theory Calculations 

To study the interfacial properties of A1 and TiN slabs 
upon approach and subsequent separation, we performed 
first-principles calculations within the framework of DFT 
employing the Vienna Ab initio Simulation Package 
(VASP)E^Jf^ VASP utilizes a plane-wave basis and pe¬ 
riodic boundary conditi ons. Projector augmented-wave 
(PAW) pseudopotential^^^^^ were used to model the po¬ 
tential between the ionic core and the valence electrons. 
Unless explicitly mentioned otherwise, the generalized 
gradient approximation (GGA) in the Perdew, Burke, 
and Ernzerhof (PBE) parametrizatioiP^ was applied to 
describe the exchange and correlation functional. Since 
GGAs often underestimate binding and adhesion ener¬ 
gies ,1^ the local-density approximation (LDA),^ which 
usually overestimates these quantities,^ was also em¬ 
ployed for comparison. Additionally, th e van der Waals 
(vdW) density functional (DE) optB86tP^^^ was used, 
which includes a nonlocal correlation term approximat¬ 
ing vdW interactions. vdW-DEs have been applied to a 
wide range of materials (e.g., see Refs. [75H83|) and have 
proven to be of good accuracy. Although vdW inter¬ 
actions should not play a major role in the investigated 
systems, the calculations are included for comparison and 
clarification. The calculation parameters were carefully 
chosen to obtain accurate total energies. An energy cut¬ 
off of 800 eV was used for the plane-wave basis. Un¬ 
less noted otherwise, the Brillouin zone sampling was 
performed using a E-centered 15 x 15 x 1 Monkhorst- 
Pack mesh.l^ Both settings allow for total energies ac¬ 
curate to 1 meV/atom. While the tetrahedron method 
with Blochl corrections^ was utilized for static calcu¬ 
lations, for relaxations a smearing of 0.11 eV using the 
first-order method of Methfessel and PaxtorP^ was se¬ 
lected. In order to relax the structures a damped molec¬ 
ular dynamics (MD) algorithm was employed, allowing 
for atomic movements until an energy convergence crite¬ 
rion of 10“^ eV was fulfilled. This damped MD scheme 
was chosen instead of the widely used quasi-Newton or 
conjugate-gradient algorithms, because these caused con¬ 
vergence problems as well as the tendency to remain 
stuck in local minima. Each converged relaxation was 
followed up by a static calculation to obtain more accu¬ 
rate total energies. Eor electronic self-consistency cycles 
a convergence criterion of 10“^ eV was used. All simula¬ 
tions were performed at OK. 


B. Simulation Model 

To model our systems we built simulation cells from a 
fee A1 slab at the bottom and a rock salt TiN slab above 
(see Eig.[^. Such cells were constructed for the low-index 
surface orientations (001), (011), and (111) of both slabs. 
Only configurations with slabs of the same surface orien¬ 


tations at the interface and without any relative rota¬ 
tions were considered. The two slabs were separated by 
a gap which is given by the vertical distance between the 
top A1 and bottom TiN layers and will be referred to as 
the “interface distance”. The vertical distance between 
the bottom A1 and top TiN layers, which is the sum of 
the interface distance and the heights of the two slabs, is 
called “slab height”. In the case of (111) slabs this height 
is measured up to the top Ti and N layer for Ti and N 
termination, respectively. Unless otherwise stated 1x1 
surface cells were used, which represent an infinitely ex¬ 
tended surface due to the periodic boundary conditions. 
The A1 slab consisted of at least seven layers, and the 
TiN slab consisted of a minimum of six Ti and six N 
layers. These thicknesses were found to be sufficient to 
converge the surface energies and to mimic bulklike fea¬ 
tures in the center of the respective slab. These system 
dimensions are in good agreement with other published 



EIG. 1: Side view of a (111) Al/TiN interface (TiN: Ti 
terminated). The simulation interface cell is indicated 
by the solid black lines. During relaxations the orange 
Al, cyan N and purple Ti atoms were kept rigid, while 
the red Al, green N and blue Ti ones were allowed to 
relax. 


The (111) TiN slab can be terminated with either Ti 
or N atoms. To investigate the st abilit y of these ter¬ 
minations a thermodynamic analysi^^^^ was performed 
to calculate the surface Gibbs free energy for the off- 
stoichiometric slabs.The surface Gibbs free energy Q 
for surface termination i without vibrational contribu¬ 
tions is given by 

^ {ELb - ( 1 ) 


where is the total energy of the slab with termi¬ 

nation i, N^- isthe number of Ti atoms in the slab, 
total energy of bulk TiN, and En is the 
total energy of a nitrogen atom. The two latter terms in 
Eq. 0 are necessary to calculate the surface energy of 
off-stoichiometric slabs. The number of off-stoichiometric 
atoms T^Ti,N defined as 
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where Nj and are the number of atoms of type j in 

the slab and in bulk, respectively. For rock-salt bulk TiN 
the fraction in Eq. ©is equal to 1. A/i^v is 

the deviation of the nitrogen chemical potential jiN from 
the molecular reference 

A^N = /^AT — (3) 

In Figure the calculated surface Gibbs free energy is 
plotted for the N- and Ti-terminated TiN (111) slabs in 
the stability range of nitrogen in TiN obtained from the 
heat of formation of bulk TihP^^ at OK, AH^ (TiN) ~ 
—3.461 eV, and the chemical potential of gas phase ni¬ 
trogen, i.e., AH^ < A/tat < 0. Fig. shows that the 
favorable termination of a (111) TiN slab depends on 
the chemical potential of nitrogen, in agreement with the 
Refs. [55] and [92l Since both cases are found in reasonable 
nitrogen concentration ranges, both terminations are in¬ 
vestigated. 

Dipole correction^ perpendicular to the interface (z 
direction) were tested for the systems but were found to 
be negligible. Atop the TiN slab a vacuum spacing of 
at least 10 A was included to decouple periodically re¬ 
peated cells in the z direction. The lattice parameters of 
the single slabs, 4.04 A and 4.254 A for A1 and TiN, re¬ 
spectively, were obtained from bulk calculations. These 
values are in very good agreement with the experimental 
lattice constants of 4.05 A and 4.265 A for A1 and TiN, 
respectively.!^ The relative error between calculated and 
experimental values is below 0.5%. For the simulation 
cells combining A1 and TiN slabs, unless otherwise stated, 
an intermediate lattice parameter of 4.144 A was used for 
the lateral xy lattice vectors to equalize the relative error 
of about 2.6% for both materials. For the z direction the 
material-specific values were kept assuming a pseudomor- 
phic interface. In reality such a combination of stretching 
and compression of thick slabs does not usually occur, 
but dislocations at the interface or an incommensurate 
contact are possible. Thus, some of the atoms on both 
sides of the interface would not be aligned perfectly, but 
rather sample slightly different local environments. For 
computational reasons, here these different local arrange¬ 
ments are assessed by considering various orientations at 
the interface as limiting cases. 

The approach of the two slabs was simulated by mov¬ 
ing the upper slab in discrete steps along the negative 
z direction and allowing for electronic and atomic relax¬ 
ations after each step. Alternatively, moving the bottom 
slab toward the upper slab or both toward each other 
would not affect the results. For the atomic relaxations 
the top TiN (three Ti and three N) and the bottom three 
A1 layers were kept fixed at bulklike distances, whereas 
the intermediate “free” ones were allowed to fully relax. 
This is depicted in Fig. for the Ti-terminated (111) 
surface orientation. For the approaching movement a 
step size of 0.2 A was used for all configurations. Before 
the slabs were brought into contact, the free layers were 
allowed to relax in order to simulate surfaces in their 



FIG. 2: Surface phase diagram for TiN. The surface 
Gibbs free energy Q [see Eq. Q] referenced to a 1x1 
surface cell of the (111) orientation is plotted vs the 
deviation Afi^ of the nitrogen chemical potential from 
its molecular reference [see Eq. ©] for N- and 
Ti-terminated (111) TiN slabs (solid lines) as well as for 
(001) and (Oil) orientations (dashed lines). 


equilibrium for the chosen lattice parameters. The sep¬ 
aration of the slabs was initiated from the equilibrium, 
i.e., the structure with the lowest energy determined dur¬ 
ing the approach. To simulate a realistic separation of 
the slabs only the topmost, rigid TiN layers were moved 
in discrete steps in the positive z direction, again allow¬ 
ing for electronic and atomic relaxations after each step. 
The choice of the step size is crucial for the separation 
process. Separation velocities allowing for an adiabatic 
behavior of the system were assumed, meaning that the 
system continuously fully adjusts during the separation 
at each step. However, this assumption should also be 
valid for velocities up to several hundred meters per sec¬ 
ond as long as these are still considerably lower than the 
material-specific speed of sound, which is above 6000 m/s 
for A1 and TiNIt is evident that the step size has to be 
small enough to mimic the adiabatic relaxation, but on 
the other hand, a smaller step size leads to increased com¬ 
putational costs. For the investigated systems a step size 
of 0.1 A was found to be a practical trade-off because cal¬ 
culations showed this value to be necessary to converge 
the final results of the simulated separation processes. 
Smaller step sizes down to 0.01 A were also considered 
for approach and separation but did not yield qualita¬ 
tively different results. Glearly, quantities such as the 
slab height corresponding to the initial material transfer 
can be determined more accurately. 

In order to study the effects of different alignments of 
the slabs at the interface, the upper slab was also later¬ 
ally placed on various sites with respect to the surface 
of the lower slab. The definitions of the configurations 
are depicted in Fig. by marking the high-symmetry 
points on the low index TiN surfaces where the next A1 
atom can be placed. In this context the interaction en¬ 
ergy Ej{z) is an important quantity, which is defined as 
the difference of the total energy of the interacting slabs 
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E(^Ai/TiN){^) at slab height 2 ; and the reference energies 
of the independent slabs, E(^ai) and 

Ei{z) = E(^Al/TiN){^) — E(AI) — E(^TiN)‘ ( 4 ) 


(001) (oil) (111) 



FIG. 3: Top view of (001), (Oil), and Ti-terminated 
(111) TiN surfaces. For each orientation the 1x1 
surface cell is presented. Filled circles indicate atoms in 
the top surface layer (Ti and N are given by large blue 
and small green circles, respectively), while empty 
circles label atoms below the top surface layer. To 
obtain a N-terminated (111) TiN surface the Ti and N 
atoms of the Ti-terminated surface have to be 
exchanged. High-symmetry points are highlighted. For 
the (Oil) TiN surface the “Ti plane” and “N plane” are 
marked by dashed lines. 


III. RESULTS AND DISCUSSION 
A. Removal of Layers from an A1 Slab 

As a first step the energy cost for removing layers from 
an A1 slab was examined for all three low-index surface 
orientations. The removal of the layers was simulated by 
placing the layers at a large distance from the slab, which 
does not allow for interactions between the slab and lay¬ 
ers. The TiN slab is not investigated here because the 
A1 slab is assumed to be mainly affected by deformations 
or material transfer within an Al/TiN interface because 
TiN forms a much more rigid lattice. The energetical 
results for the removal of the top A1 layer are given in 
Table [I| These removal energies are calculated for simula¬ 
tion cells using the bulk lateral lattice parameters as well 
as the modified ones used for the Al/TiN simulation cell. 
For the modified lattice parameters the removal energies 
are typically overestimated by about 5%-10%, meaning 
that it is actually easier to remove layers from the equi¬ 
librium structure. The removal energy for the modified 
A1 slab is increased because the lateral stretching causes 
a vertical compression of the slab if relaxations are al¬ 
lowed. This compression occurs to minimize the volume 
change and locally strengthens the bonding of the surface 
layers. This effect is strongest for the top surface layer, 
which moves about 0.24 A towards the rigid part and be¬ 
comes weaker for the subsurface layers; for example, the 
fourth layer is only shifted by about 0.08 A. 


The influence of compressive and tensile stress on the 
removal energies of the top A1 layer is illustrated in Fig.[^ 
for the three low-index surface orientations. The data 
points for the (001) and (111) surfaces follow a similar 
trend, whereas the behavior of the (Oil) surface clearly 
deviates. This difference occurs probably due to the 
openness of the (Oil) surface and the significant im¬ 
pact of relaxations. The influence of stress, found f or all 
surfaces, supports the notion of stress-assisted wear,C^E^ 
which states the possibility of a reduction of the activa¬ 
tion barriers for the detachment of atoms from a struc¬ 
ture due to stress. Furthermore, different approxima¬ 
tions for the exchange correlation functional were tested. 
As expected it was found that LDA and the vdW-DF 
optB86b yield larger removal energies by about 15%- 
20%, where LDA typically gives values larger by a few 
percent than the vdW-DF (see Table |I|. 


TABLE I: Energy costs to remove the top layer from 
an A1 slab for the (001), (011), and (111) surface 
orientations using PBE, LDA, and vdW-DE optB86b. 
The removal energies are given in eV per 1x1 surface 
cell. aAi is the A1 bulk lattice parameter, whereas 
^AijTiN corresponds to the modified Al/TiN interface. 



(001) 

(Oil) 

(111) 

PBE (aAi) 

1.08 

1.78 

0.78 

PBE (aAi/TiN) 

1.16 

1.79 

0.87 

LDA (aAi/TiN) 

1.33 

2.02 

1.02 

vdW (aAi/TiN) 

1.27 

1.89 

1.00 



EIG. 4: PBE energy costs to remove the top A1 layer 
for the (001), (011), and (111) surface orientations. The 
lateral effects of stretching and compression of the 1x1 
surface cell on the removal energies are shown. The A1 
bulk lattice constant is used as a reference value at 0%. 

The vertical line indicates the intermediate Al/TiN 
interface lattice parameter, while the other solid lines 
are given to guide the eye. 
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B. Lateral Alignments at the Al/TiN Interface 

Effects of various lateral alignments of the slabs at the 
interface (see Fig. were investigated for the different 
surface orientations. These studies revealed the strong 
dependence of equilibrium properties such as adhesion 
energies and the equilibrium distances on the chosen con¬ 
figuration. The calculated interaction energies [Eq. @] 
of relaxed interfaces are shown in Figs. |^a)-(d) for slab 
heights around the energy minima, which are equivalent 
to the adhesion energies for each alignment. In general, 
the top A1 atoms prefer the proximity of N atoms over 
Ti. The bonding situation will be discussed in more de¬ 
tail in the following paragraphs. From an energetical 
point of view material transfer between the slabs should 
be possible only if the energy cost to remove layers is 
compensated for. Thus, the energy gain due to adhesion 
has to be larger than the energy cost to remove one or 
more layers. This argument is sketched in Figs. ia)- 
(d) by including a horizontal line at the negative value 
of the A1 removal energy for each surface orientation. It 
has been observed experimentally that metal-ceramic in¬ 
terfaces with weak and strong interfacial adhesion break 
upon stres s at the interface and in bulk areas, respec¬ 
tively We find that the four surfaces investigated 
exhibit essentially different behavior. The adhesion en¬ 
ergies and the equilibrium distances, i.e., the interface 
distances at the minimum of each energy curve, depend 
strongly on the surface orientation as well as on the align¬ 
ment at the interface. In the case of the (111) surfaces all 
configurations should lead to the removal of at least one 
A1 layer. For the (Oil) surfaces this is the case only for 
three alignments (see Fig. |^, Al/N (top), Al/TiN (hol¬ 
low), and Al/N (bridge). In contrast, for the (001) sur¬ 
faces no material transfer should occur since for all cases 
studied the energy to remove one A1 layer is larger than 
the adhesion energy. 

As mentioned above, in reality, surfaces with a bulk 
lattice mismatch are usually not perfectly aligned at an 
interface. Consequently, not all atoms are placed on the 
same contact site; therefore, the interfacial properties 
such as the adhesion energy are an average of the actually 
occupied sites. The configurations presented here, how¬ 
ever, constitute limiting cases of perfectly aligned sys¬ 
tems, such that the properties of real interfaces should 
be found within these boundaries. 

Generally, relaxation effects have to be accounted for 
to obtain the correct equilibrium values of the adhesion 
energy and the interface distance as well as to predict 
the occurrence of material transfer. A comparison be¬ 
tween the relaxed and static results is given in Fig. for 
the (111) surfaces. For rather closed TiN surfaces, such 
as the (001) and Ti-terminated (111) orientations [see 
Fig.§ a)], relaxations typically cause only small changes 
in the equilibrium quantities of the interface. Hence, 
computationally “cheap” static calculations give good es¬ 
timates, unless pronounced changes in the structure of 
the A1 slab occur. This is, for example, the case for 


the Al/Ti (hollow) alignment of the (111) Al/TiN (Ti- 
terminated) interface, since the inter facial A1 atom re¬ 
laxes towards the energetically more favorable fee con¬ 
tact site. In the case of the more open (Oil) surface, 
relaxations show more pronounced effects for all align¬ 
ments and should be taken into account. Nevertheless, 
the energy hierarchy and the prediction of the occur¬ 
rence of material transfer are not affected for all align¬ 
ments with the exception of the Al/TiN (hollow) case. 
Again, the Al/TiN (hollow) interface behaves differently 
because the relaxed structure of the A1 slab is modi¬ 
fied by the approaching TiN slab. In more detail the 
interfacial A1 layer is moved to the Al/N (bridge) site, 
which is the most favorable alignment. This movement 
of about 0.8 A occurs mainly in the lateral plane. The 
free subinterface layers are shifted to gradually compen¬ 
sate the change in the stacking between the fixed layers 
at the bottom of the slab and the interfacial layer. These 
shifts range approximately between 0.2 and 0.6 A. For 
the cases discussed so far, except for the hollow align¬ 
ments, relaxations showed rather small effects on the 
equilibrium quantities. In contrast, all alignments of the 
(111) Al/TiN (N-terminated) interface are crucially af¬ 
fected by relaxations [see Fig. [^b)]. The adhesion ener¬ 
gies are strongly increased, and the energetical hierarchy 
of the alignments is altered. Furthermore, while static 
calculations suggest the absence of material transfer, re¬ 
laxations predict its occurrence for all tested alignments. 

For a better understanding of the energetically pre¬ 
ferred configurations at the interface, layer-projected 
densities of state (DOSs) and differences in charge den¬ 
sities are examined. Layer-projected DOSs are displayed 
in Fig. for the two alignments Al/N (bridge) and 
Al/Ti (top) as well as the isolated slabs of the (Oil) sur¬ 
face orientation. This surface orientation has been chosen 
because it exhibits a large spread in adhesion energies for 
different alignments. Additionally, the occurrence of ma¬ 
terial transfer should depend on the alignment. In Fig. 
“interface (surface) layers” indicate the first layers of Al, 
Ti, and N immediately at the interface (surface), whereas 
“subinterface (subsurface) layers” mean the next layers of 
Al, Ti, and N moving deeper into both materials. Further 
layers are not presented because they exhibit only minor 
differences with respect to the subinterface layers. The 
DOSs of the shown alignments display distinct features. 
For the Al/Ti (top) case, where Ti is the next interfacial 
neighbor of the top Al atom, the Al DOS is almost not 
affected by the interface. Only a small accumulation of 
sp states just below the Fermi energy and a depletion of s 
states at the edges of the DOS are found for the interface 
layers with respect to the other layers. The N sp states 
are shifted closer to the Fermi energy for the interfacial 
layer, and in particular, the Ti d states exhibit more oc¬ 
cupied states at the Fermi energy. These changes indicate 
a weakly covalent bonding between the Al sp states and 
the Ti d states. Furthermore, the DOS is very similar 
to the case of the isolated Al and TiN slabs. This also 
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FIG. 5: Calculated PBE interaction energies of the relaxed Al/TiN interface for the (111) Ti-terminated, (111) 
N-terminated, (Oil), and (001) surface orientations. Various lateral alignments of the two slabs are considered (see 
Fig.§. The horizontal, dashed orange lines give the energy costs to remove at least one layer from an A1 slab of the 

corresponding surface orientation. 


shows the weak interaction for the Al/Ti (top) interface. 
On the other hand, for the Al/N (bridge) configuration, 
where the uppermost A1 atoms are closer to N across the 
interface, the A1 DOS is changed in a more pronounced 
way. The sp states in the interface layers are partially 
shifted to lower energies, resulting in a pronounced peak 
at about -8 eV and a few minor ones around -7 eV. The 
N p states around -5 eV are broadened in the interfacial 
layer, resulting in common peaks with A1 states roughly 
between -6 eV and -8 eV. These effects at the interface in¬ 
dicate a hybridization of A1 and N sp states and explain 
the stronger adhesion due to covalent interaction. The 
interfacial Ti states are only slightly affected, exhibiting 
a few more occupied states at the Fermi level. 

In addition to the DOS, charge densities at the inter¬ 
faces are investigated and presented for the same align¬ 
ments of the (Oil) surface. For a better visualization the 
differences of charge densities Pdiff between the Al/TiN 


interface and the isolated, independent A1 and TiN slabs 
are presented in Fig. The charge-density difference 
Pdiff is defined as 

Pdiff = PAllTiN — {pAl + pTiN), (5) 

where pAijTiN is the charge density of the interface, while 
Pai and pTiN represent the charge densities of the iso¬ 
lated slabs. Both displayed alignments result in a rather 
continuous charge accumulation between A1 and Ti at the 
interface, suggesting a bonding [see Figs. a) and[^c)]. 
For the Al/N (bridge) configuration an additional charge 
buildup occurs between the interfacial A1 and N atoms, 
which indicates covalent contributions to the bonding due 
to the more localized and directional character of the ac¬ 
cumulation [see Fig. id)]. These findings support the 
DOS arguments made in the previous paragraph. 
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(a) (111) Ti-terminated 


(b) (111) N-terminated 


FIG. 6: Calculated PBE interaction energies of the Al/TiN interface for the (111) Ti- and N-terminated surfaces. 
Solid and dashed lines indicate results of relaxed and static calculations, respectively. The slab heights on the lower 
X axis are valid for static and relaxed calculations, whereas the interface distances on the upper x axis refer only to 
the static calculations. Various lateral alignments of the two slabs are considered (see Fig. |3). The horizontal solid 
orange lines give the energy costs to remove at least one layer from an A1 slab. 


C. Approach and Separation of A1 and TiN Slabs 

The energetical argument on material transfer pre¬ 
sented above can be tested by “slowly”, i.e., using small 
discrete steps, approaching and subsequently separat¬ 
ing the slabs. The energetical results of such loops are 
depicted in Figs, [^a)-(d) for different configurations; 
the respective energies are presented in Table |TI[ The 
green curves with their data points indicated by pluses in 
Figs. I^a) -(d) give static potential-energy curves, where 
all atoms were kept rigid for each selected interface dis¬ 
tance. For large interface distances this curve shows the 
limiting case of separated, independent slabs. For ever- 
shorter distances the effect of relaxation becomes im¬ 
portant, and the actual energies deviate from the green 
curves. The blue curves with their data points marked 
by crosses show the interaction energies of the approach¬ 
ing slabs including atomic relaxations after each discrete 
step. For some of these cases we find rather large jumps 
which are either due to material transfer between the 
slabs, namely, from A1 to TiN, or due to the A1 slab ex¬ 
panding into the space between the slabs. Finally, the red 
curves with their data points displayed by circles indicate 
the interaction energies of the subsequent separation of 
the slabs, again including atomic relaxations after each 
step. These curves are also not completely smooth but 
display some kinks or smaller jumps mainly due to the 
breaking apart of the Al/TiN slab into two separated 
ones. When material transfer takes place, these curves, 
of course, do not approach the green ones, even at large 
slab separations. 

For the Ti-terminated (111) surface orientation 
potential-energy curves are presented in Fig. [^a) for the 


two extremal alignments, Al/Ti (hep) and Al/Ti (top), 
which show the highest and lowest adhesion energies. 
As expected from the energetics, material transfer oc¬ 
curs during separation, and both systems end up in an 
energetically more favorable configuration compared to 
the initial setup. In particular, one and two A1 layer(s) 
for Al/Ti (top) and Al/Ti (hep), respectively, are trans¬ 
ferred. This discrepancy in the number of transferred 
layers cannot be explained from the energetics but could 
stem from the different equilibrium interface distances. 
Compared to that of the hep alignment, this distance is 
significantly increased by almost 20% for the top config¬ 
uration, hindering the interaction between TiN and the 
subinterface A1 layer. For the Al/Ti (hep) configuration 
snapshots of the structures during approach and sepa¬ 
ration are presented in Fig. [OT During the approach 
at a slab height of about 33.6 A a large drop in interac¬ 
tion energy occurs for the Al/Ti (hep) alignment [see 
Fig. [^a)] due to material transfer of the topmost A1 
layer to the TiN slab [see Fig. [^b)]. This is not the 
ground state since a transfer of two layers would yield an 
even lower total energy. At this distance, the transfer of 
the second A1 layer is hindered by an energy barrier of 
about E })2 ~ 324 meV, which is significantly larger than 
for the first layer alone, « 122 meV. Upon further 
approaching, at a slab height of about 32.8 A, a slight 
kink occurs [see Fig. |^a)] because the A1 slab is ex¬ 
panded into the space between the slabs [see Fig, [l^c)]. 
For a further approach, the interaction energy follows 
an essentially parabolic curve until the minimum energy 
is reached [see Figs, [^a) and [T^d)]. The subsequent 
separation is started from the equilibrium structure at 
the interaction energy minimum. At first the red in¬ 
teraction energy curve lies on top of the blue one [see 
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FIG. 7: Layer-projected DOSs from PBE calculations of the isolated A1 and TiN slabs as well as of the (Oil) 
Al/TiN interface for the Al/Ti (top) and Al/N (bridge) alignments. The Fermi energy is shifted to OeV. 


Fig. [^a)], meaning that the A1 slab becomes extended 
again [see Fig. |l^f)]. At a slab height of about 33.1 A 
the two curves for approach and separation start to devi¬ 
ate [see Fig.j^a)] when the Al/TiN compound separates 
[see Fig. 10 [g/]. Two A1 layers stick to the TiN slab and 
form a stable configuration. This behavior during the 
complete loop is typical for all cases exhibiting material 
transfer. While the A1 slab is strongly affected by the 
approach of the TiN slab, almost no changes in the TiN 
structure are observed. The more pronounced impact on 
the A1 slab is not surprising when considering that TiN 


forms a much more rigid lattice than Al. This claim is not 
entirely valid for the N-terminated (111) TiN slab, which 
will be discussed in the following paragraph. Using the 
finally stable state (TiN plus two Al layers) as a starting 
configuration for a new loop of approach and separation 
versus an Al slab yields a reversible cycle. This should be 
kept in mind when one is interpreting, for example, AFM 
experiments. Upon the first contact between the tip and 
a particular spot on a surface material transfer might oc¬ 
cur, which in turn changes the contact properties and 
forces between the tip and the surface. However, further 
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TABLE II: Equilibrium interface distances, adhesion energies, energy costs to remove layers from the A1 slab, and 
number of transferred A1 layers for various interface configurations. Eor the (111) orientation Al/Ti and Al/N 

denote the Ti- and N-terminated surfaces, respectively. 


(001) Al/N (top) 

Equilibrium 
interface distance [A] 

2.06 

Adhesion energy 
[eV/interface cell] 
-0.61 

Removal energies 
[eV/interface cell] 

1.16 

Material transfer 
[Al layers] 

0 

(Oil) Al/N (bridge) 

1.39 

-2.09 

1.35 

2 

(Oil) Al/Ti (top) 

2.77 

-0.73 

1.35 

0 

(111) Al/Ti (hep) 

2.22 

-1.78 

0.80 

2 

(111) Al/Ti (top) 

2.67 

-0.94 

0.80 

1 

(111) Al/N (hep) 

1.04 

-1.90 

0.80 

2 

(111) Al/N (top) 

1.87 

-1.38 

0.80 

1 
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(b) 

Al/Ti 

(top): 

N-plane 


(c) 

Al/N 

(bridge): 

Ti-plane 


(d) 

Al/N 
(bridge): 
N-plane 


EIG. 8: Charge-density differences pdiff [see Eq. 0] of 
the (Oil) Al/TiN interface, pdiff was obtained from 
PBE calculations for the relaxed equilibrium 
configurations of (a) and (b) the Al/Ti (top) alignment 
and (c) and (d) the Al/N (bridge) alignment. The 
charge-density difference of each alignment is plotted 
for the Ti plane and the N plane (recall Eig. ^ for 
values from -0.2 (solid blue, deficit) to 0.2 (solid red, 
accumulation) electrons/A^. Color code: Al, orange; Ti, 
violet; N, cyan. 


encounters on the same spot should then be within the 
reversible cycle and lead to the same response. 

The N-terminated (111) orientation is, in some re¬ 
spects, very similar to the Ti-terminated one. As pre¬ 
dicted, both configurations yield material transfer for all 
tested alignments (see Eig. |^. However, as explained 
above, static calculations completely fail to describe the 
equilibrium quantities of the N-terminated case, whereas 
they result in good estimates for the other orientations. 
This discrepancy is due to the behavior of the interfacial 
N layer for N-terminated (111) TiN. In the absence of 


the counter Al slab, the surface N layer is closely bound 
to the next Ti layer at a distance of 0.84 A, while in 
contact with an Al slab the distance grows to 1.47 A at 
the equilibrium configuration. This behavior is crucial 
for the energetics and can be captured only when relax¬ 
ations are included. The inter facial N layer is actually 
closer to the next Al layer with a distance of 1.11 A than 
to the next Ti one. Due to this result the possibility 
of a diffusion of the interfacial N layer into the Al slab 
was investigated. Eor all alignments with the exception 
of Al/N (top), no energetically favorable configurations 
were found. However, for the Al/N (top) alignment the 
exchange of the interfacial Al and N layers and a sub¬ 
sequent relaxation of the system yield a favorable state 
by about 683 meV, which is also about 235 meV lower 
than the previously found minimum for the Al/N (hep) 
alignment. In this favorable configuration an Al-N-Al 
trilayer is formed, showing the wurtzite structure, which 
is typically observed in aluminum nitride crystals. Erom 
the thermodynamical point of view diffusion seems to be 
possible. Of course, for the full picture reaction paths 
and energy barriers have to be considered. 

The (Oil) surface orientation also presents an inter¬ 
esting case because due to the energetic results [see 
Eig.j^c)], material transfer is expected for only the align¬ 
ments Al/N (bridge), Al/N (top), and Al/TiN (hollow). 
As an example, one can see from the loops given in 
Eig.j^c) that Al/N (bridge) shows a favorable configura¬ 
tion after separation corresponding to the transfer of two 
Al layers, whereas the Al/Ti (top) case is reversible upon 
approach and separation without any material transfer. 

Einally, for the (001) surface orientation material 
transfer is not expected for any of the alignments. Among 
all cases Al/N (top) has the largest adhesion energy; 
therefore, if a material transfer occurs, it will happen 
for this case. However, since the energy cost for the re¬ 
moval of an Al layer exceeds the adhesion energy, no 
material transfer is observed [Eig. [^d)]. The deviation 
of the curves for approach and separation around 26.5 A 
slab height occurs due to the expansion of the Al slab 
upon separation until the interface breaks apart and re¬ 
laxes into the initial Al and TiN slabs. 

In the literature some publications on tensile test sim- 
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FIG. 9: Calculated PBE interaction energies [see Eq. for the approach and subsequent separation of A1 and TiN 
slabs for (111) Ti-terminated, (111) N-terminated, (Oil), and (001) surface orientations. The alignments follow the 

definitions in Eig. 


ulations of Al/TiN interfaces can be found, where the 
separation is achieved by increasing the size of the whole 
simulation cell in one direction in discrete steps includ¬ 
ing interim relaxations. Liu et and Zhang et al.^^ 
investigated the Al/TiN (111) and (001) interfaces, re¬ 
spectively. Liu et al. obtained similar results with re¬ 
spect to material transfer for the hep alignment at the 
(111) surface for both terminations but did not examine 
any further alignments of Al and TiN slabs at the inter¬ 
face. Zhang et al. studied the Al/N (top) configuration 
of the (001) interface and, in contrast to our work, found 
a material transfer of the top Al layer. This discrepancy 
could stem from the different simulation approaches and 
computational details. However, we repeated these cal¬ 
culations using a setup for the separation of the slabs 
similar to that of Zhang et al. and did not find any ma¬ 
terial transfer. Additional simulations for the different 
setups used by Zhang et al. and in the present investi¬ 
gation testing the influence of varying step sizes also did 
not lead to a material transfer. 


D. Comparison of Surface Energies 


The behavior of the different surface orientations can 
also be discussed from the surface energy’s point of view. 
The sur face energies for Al and TiN slabs are presented 
in Table III It has to be noted that for (111) TiN the sur¬ 
face energy depends on the termination and the chemical 
potential of nitrogen. Here the lowest possible value for 
the surface energy is used, which is achieved by the N- 
terminated surface at Ajun = 0 (see Eig. |^. The value 
of the surface energy for the Ti-terminated surface at 
A/i 7 v = 0 is about three times larger, but its minimum 
is comparable to that of the N-terminated case. Eor the 
(001) and (011) orientations the surface energy is inde¬ 
pendent of the chemical potential.!^ As shown in Ta¬ 
ble |HI[ the Al surfaces always exhibit a smaller surface 
energy than the TiN ones. The differences between Al 
and TiN are pronounced for the (011) and (111) orienta¬ 
tions. Eor (001), however, the surface energies are rather 
comparable. Material transfer can be seen as a measure 
of surface-energy minimization by creating a new ener¬ 
getically “cheap” surface and covering an “expensive” 
one with it. This argument provides a hint about which 
surface orientations may favor material transfer. How- 
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FIG. 10: (a)-(d) Approach and (e)-(h) separation of 
(111) A1 and (111) TiN slabs. Al, Ti and N are colored 
in red, blue, and green, respectively, (d) and (e) show 
the structure at the relaxed equilibrium distance. The 
steps are defined by the slab height. 


ever, for the full picture also other contributions such as 
the interaction energy, which is influenced additionally 
by the alignment of the slabs, also have to be considered. 
For example, the surface-energy argument would suggest 
the possibility of material transfer for (001) and cannot 
explain why only some (Oil) configurations exhibit this 
feature. 


E. Assessment of Computed Results 

To validate the information presented above various 
additional tests were performed. The results will be pre¬ 
sented for the Ti-terminated (111) Al/Ti (hep) configu¬ 
ration. First of all, finite-size effects are a major concern. 
Thus, the size of the simulation cell was increased later¬ 
ally up to a 3x3 surface cell and vertically up to a 19- 
layer Al slab. Also, intermediate Al slab thicknesses were 
examined. The TiN slab was not extended vertically be¬ 
cause it is almost not affected by the approach of the Al 
slab. In the case of laterally magnified simulation cells 
the number of k points was decreased accordingly, e.g., 
for a 3x3 surface cell, a 5x5x1 mesh was used. For all 
tested systems the equilibrium interface distances, adhe¬ 
sion energies, and energy costs to remove Al layers were 
found within about 2% of the values given above. The 
energies are referenced to 1x1 surface cells. Particularly, 
the results on material transfer were not affected, mean¬ 
ing that the number of transferred Al layers was not al¬ 
tered. For the 3x3 surface cell the effect of fluctuations 
at the surface was tested by moving one atom out of the 
surface plane at several interface distances before mate¬ 
rial transfer occurs. These tests resulted in the transfer 
of entire layers too since the shifted atom either relaxed 
back into its originating slab or was transferred together 
with the rest of the layer. 

Furthermore, the effect of the chosen lattice parame¬ 
ters was investigated. The simulations were repeated us¬ 
ing the lattice constants of pure Al and TiN for the lateral 
lattice parameters of the simulation cell. The equilib¬ 
rium interface distance again changed by only about 2%. 
Although the adhesion energies were altered by about 
4%, the removal energies were affected in a similar way, 
resulting in the same material transfer. Moreover, the 
influence of other approximations for the exchange cor¬ 
relation functional was tested as already discussed for 
the removal energies above. The results are presented 
in Fig. EH The adhesion energies were enhanced sim¬ 
ilar to the removal energies, again producing the same 
results for material transfer. Using the vdW functional, 
the interaction between the slabs started at larger inter¬ 
face distances. This behavior is expected because of the 
nonlocal correction added in the vdW functional. 


TABLE III: Surface energies (in eV/A^) of the Al and 
TiN slabs for the (001), (011), and (111) surface 
orientations. In the case of the (111) TiN surface the 
N-terminated one at A/j^n = 0 is given here because it 
exhibits the lowest surface energy of all (111) TiN 
surfaces (see Fig. [^. 


(001) (Oil) (111) 

Al 0.058 0.064 0.052 

TiN 0.087 0.174 0.094 


IV. CONCLUSION AND OUTLOOK 

Al/TiN interfaces were examined in detail by inves¬ 
tigating the contact between Al and TiN slabs showing 
low-index surface orientations within the framework of 
density functional theory. Moreover, these contacts were 
established for various lateral alignments of the slabs at 
the interface. It was shown that interfacial properties 
such as the adhesion energy and the equilibrium struc¬ 
ture sensitively depend on the given configuration. This 
behavior can be qualitatively explained by comparing the 
densities of state and the charge densities of different con- 
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FIG. 11: Comparison of calculated interaction energies 
for the Ti-terminated (111) Al/Ti (hep) configuration 
and various exchange correlation functionals, namely, 
PBE, LDA, and opbB86b. The relaxed energies 
represent the separation of the slabs starting from the 
equilibrium configuration. The differences between 
static and relaxed curves at large heights occur due to 
material transfer. 


given alignment at the interface has to be considered as 
well. All results were tested for various computational 
setups such as different sizes of the investigated system 
or several approximations for the exchange correlation 
functional. While properties such as the removal and ad¬ 
hesion energies depend on these settings to some degree, 
the results for material transfer are not affected. 

The method used in this work can be, in principle, ap¬ 
plied to any pair of materials. However, complex mate¬ 
rials or pairs with an unfavorable bulk lattice mismatch 
may need very large simulation cells to be considered, 
which means high computational demands, in order to 
preserve the translational symmetry and to keep the dis¬ 
tortions at an acceptable level. Furthermore, larger cells 
also allow the inclusion of additional features. For exam¬ 
ple, the distortions due to the lattice mismatch can be 
minimized, dislocations as well as quasi-incommensurate 
contacts can be modeled, and even roughness could be 
included to some degree, e.g., by using stepped surfaces 
or a regular grid of small asperities. 


figurations because distinct bond situations are revealed 
at the interface. Furthermore, the approach and subse¬ 
quent separation of A1 and TiN slabs were simulated to 
study the effect on the slabs, especially the possibility of 
material transfer. The transfer of material from an A1 
toward a TiN slab was observed for interfacial configura¬ 
tions, which exhibited a larger adhesion energy than the 
energy cost to remove layers from the A1 slab. This is in 
agreement with the observation that metal-ceramic in¬ 
terfaces break at the interfa ce or in bulk areas according 
to their interfacial adhesion.The removal energy for 
A1 layers was found to depend on tensile or compressive 
stress. In all systems showing material transfer one or 
two layers of A1 stick to the TiN slab after the separa¬ 
tion and form an energetically favorable compound with 
respect to the initial configuration. The differences in 
surface energies between the slabs are not sufficient to 
explain the occurrence of material transfer because the 
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